Symmetry Breaking in Stock Demand 
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Scale-free distributions and correlation functions found in financial data are reminiscent of the 
scale invariance of physical observables in the vicinity of a critical point. Here, we present empirical 
evidence for a transition phenomenon, accompanied by a symmetry breaking, in the investors' 
demand for stocks. We study the volume imbalance Q. — difference between the number of shares 
traded in buyer-initiated and seller-initiated trades in a time interval At — conditioned on E which 
is defined as the local first moment of Q. in At. We find that the conditional distribution P(r2|E) 
undergoes a qualitative change in behavior as E increases beyond a critical threshold Ec. For 
E < Ec, P(S1|E) displays a maximum &t Q, = 0, i.e., trades in At are equally likely to be buyer 
initiated or seller initiated. For E > Ec, f2 = becomes a local minimum and two new maxima 
Q.+ and n_ appear at non-zero values of Q., i.e., trades in At are either predominantly buyer 
initiated or predominantly seller initiated. We interpret these results using a Langevin equation 
with multiplicative noise. 
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Phase transitions of systems from an 'ordered' phase 
into a 'disordered' phase are closely linked to symmetry 
breaking. For example, in an Ising ferromagnet above 
its critical temperature, the most probable state (zero 
net magnetization) possesses the symmetry that leaves it 
invariant under the flipping of each spin. For tempera- 
tures below the critical value, this symmetry is broken, 
and there is a preferred direction. 

Here, we present evidence for an analogous transition 
phenomenon in a financial context [0-j7j. Specifically, 
we study the statistical properties of investors' demand 
for stocks — quantified as the imbalance in the num- 
ber of shares transacted by buyers and sellers over a 
time interval At. We analyze the probability distribu- 
tion of demand, conditioned on its local "noise" inten- 
sity (a variance- like parameter E defined below) . We find 
that for intensities smaller than a critical value Sc, the 
most probable value of demand is approximately zero — 
neither buying nor selling behavior dominates. For in- 
tensities larger than that critical value (E > Ec), two 
most probable values emerge that are symmetric around 
zero demand, corresponding to two distinct "phases" — 
excess demand and excess supply Under such con- 
ditions, the market behavior is either mainly-buying or 
mainly-selling, spending almost equal amount of time in 
each state. In other words, exchanging every "buy" with 
a "sell" gives the same state below the critical noise in- 
tensity, whereas above this threshold, the symmetry of 
this exchange is broken. 

In classic critical phenomena, the qualitative change in 
behavior accompanying a phase transition can be formal- 
ized in terms of the extrema of a phenomenological po- 
tential, or equivalently in terms of the cxtrcma of the cor- 
responding probability distributions [Hjoj . We first follow 
the latter approach and study the behavior of the prob- 
ability distribution of demand. Using transactions and 



quotes data for the 116 most-actively traded stocks pO[ , 
we quantify the demand for a stock by calculating the 
"volume imbalance" over a time interval At, defined to be 
the difference between Qb, the number of shares traded 
in buyer-initiated trades, and Qsj the number of shares 
traded in seller-initiated trades in At |lll-[l^. 
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Here, the indicator = 1 for buyer-initiated trades 
(buy trades) and = — I for seller-initiated trades (sell 
trades) ||l7|] , qi is the number of shares traded in trans- 
action i, and N = Nj\t denotes the number of trades in 
At 0,0. 

Our analysis of the (unconditional) probability distri- 
bution P(n) for each stock shows a single peak around 
f2 = 0. Since previous work shows that the distribution 
of Qi has divergent variance |]l9| , ^ , we quantify the noise 
intensity by computing the 'local deviation' pl| ], defined 
as the centered first moment. 



E(t) = {\qiai - {qiai)\) , 



(2) 



where (...) denotes 'local' expectation values computed 
from all trades in the time interval At. Next, we examine 
the behavior of the conditional distribution P(ri |E) of 
for a given value of the local deviation E for At = 15 
min, [Fig. 1(a)]. For small E, we find that P(f7|E) is 
single peaked displaying a maximum at Q = 0. When E 
exceeds a critical threshold Ec, the behavior of P(f7|E) 
undergoes a qualitative change, and is double peaked with 
two new maxima appearing at non-zero values, fl+ and 
ri_, symmetric around zero. Figure 1(a) also shows that 
the separation between the two maxima increases with 
E. 
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This quaUtative change in the behavior of P{^) im- 
plies that for S < Sc, the most-probable value of demand 
is approximately zero, and possesses the symmetry that 
leaves the most probable value invariant under the oper- 
ation — > — f2, or at the microscopic (trade) level, under 
the operation B ^ S of changing every buyer-initiated 
trade B, to a seller- initiated trade S. For E > Ec, the two 
most probable values fl± are non-zero, and the Q — O 
(B S) symmetry is broken. In other words, while for 
E < Ec buy and sell trades are equally probable in each 
time interval (zero demand), for E > S^, trades in each 
time interval are either mostly buy trades (excess de- 
mand) or mostly sell trades (excess supply) giving rise to 
non-zero values of . Identical results can be obtained 
by conditioning P{^1) on the total trade volume in At, 

g(t) = Qb + Qs. 

Our finding is analogous to phase transition phenom- 
ena in physical systems, where the behavior of the system 
undergoes a qualitative change at a critical threshold of a 
control parameter T. In such systems, the change in be- 
havior can be quantified by an order parameter which 
is identically zero for values of T below (or above as the 
case may be) a certain critical value T^, and becomes 
nonzero as T crosses Tc. In our problem, the "order 
parameter" ^P, can be identified by the location of the 
maxima il± of P{^)- Figure 1(c) shows that the change 
in \1/ as a function of E is described by 
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[E < Ec 
fE > E, 



(3) 



with /? « 1 g]. 

In the mean-field theory of critical phenomena 
(Landau-Ginzburg theory) , the qualitative change in be- 
havior of the system is attributed to the changes in sym- 
metry of the underlying potential |^,^. In the follow- 
ing, we pursue an analogous approach to understand 
our empirical results. Since the transition behavior that 
we find occurs with change in noise intensity, we follow 
an approach similar to those used to understand non- 
equilibrium phase transition phenomena p3|-p5t. We 
start with expressing the dynamics of Q through a de- 
terministic differential equation. 



dn = hx{n)dt, 



(4a) 



where A is a parameter quantifying the coupling of the 
system to its environment. Letting A fluctuate randomly 
with noise intensity cr, Eq.(|4^) becomes, in general, a 
stochastic differential equation with multiplicative noise. 
Since the form of h\(ri) is not known, and since it is 
not a priori clear if A is an observable, we describe the 
dynamics of through ||23] 



dfl = u{n)dt + av{n)dWt 



(4b) 



where u = u{fl) is the drift term, v = v(Q) reflects 
the effects of multiplicative noise, dWt is the standard 
Wiener differential satisfying (dWtdWt') = S{t - t')dt, 



and t he par ameter cr quantifies the intensity of the noise 
term |26|j2^ . The functions u and v may estimated from 
the data as a "local" mean and standard deviation of 
dO. l2q|. The Fokker-Planck equation corresponding to 



Eq. (4b) which describes the evolution of the transition 
probability density 11 = P{Vl,t\VL{t — 0)) is 
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The stationary probability density P{^) from Eq. 



^(^)-^exp(A r^Ld. 



(5) 
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where A is the normalization constant, and P is assumed 
to be normalizable. Equation (H) can be rewritten as 



P{n) = A exp 



(J) 



where the function V{^) takes the meaning of a effective 
"stochastic" potential given by [p3| 



v{n) 



-irdx — a \nv 
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Thus, the extrema of the probability density can be 
translated into the extrema of the underlying stochastic 
potential. 

As the noise intensity of the environment (a) changes, 
the potential could change shape, acquire new min- 
ima, and consequently the system could undergo drastic 
changes in behavior p9| . Can the transition behavior we 
find empirically for PjTl\J^) be understood in this frame- 
work? To address this question, we must examine the 
shape of the potential V{^1) for different values of E, 
which monotonically depends on cr |30|. 

In order to study V(fl) for different cr empirically, we 
first extract the functions u{Q) and v{fl) and analyze 
their behavior for different values of E. From Eq. (^), 
u is the drift term, which for a given E, can be esti- 
mated by computing the equal-time expectation value of 
the change Af2 = n{t + At) — il{t) for a given fl, 



u « (Arj)(o,s) 



(9) 



Similarly, the product av can be estimated from the "lo- 
cal" deviation 



av 



(|Ar!- (Al])(o,s)|)(oj 



(10) 



Figures ^( a) and (b) show u{^l) and v{Vl) for three dif- 
ferent values of E. Clearly, the functional form u{Q) does 
not vary with E and is consistent with a linear behavior 
for all E. Figure |^(b) shows that, for small E, (jv{Vl) 
is approximately flat, whereas for large E, crw acquires 
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a marked 'peak' around = 0. Except for the small- 
est S, the functional forms av{fl) for different E seem 
to be consistent within a multiplicative factor (related to 
a) H. 

Next, we shall analyze the extrema of the stochas- 
tic potential V{fl) for different values of E. From 
Eqs. @j and (^), the extrema of correspond to the 

roots of the function 



Fin) 



dvjn) 



dv 
dh' 



(11) 



Figure||(a) shows F{Q) for three different S. For E < Ec, 
the function F{n,) displays only one root at = 0. Near 
E ~ Ec, F{n) displays an inflexion point at — 0. For 
E > Ec, the one root branches into three: = 0, 17+, 
and r2_. 

Integrating the stochastic "force" F{n), we next find 
the potential V{Q). Figure ^ shows the behavior of 1^(17) 
with different E. For E < Ec, the potential has only one 
minimum at — 0, which is consistent with one maxi- 
mum for P(ri) which we find. For E « Ec, the poten- 
tial appears almost flat: the existing minimum begins to 
change into a (unstable) local maximum and displays an 
inflexion point at E = 0. For S > Ec, the potential dis- 
plays two clear minima at flj- and a local maximum at 
n = 0, consistent with the bimodal nature of the distri- 
bution P{n) for E > Ec. 

In summary, we investigate the dynamics of the de- 
mand 17 by examining the distribution of volume imbal- 
ance O for changing market conditions quantified by the 
local deviation E. We find that the distribution P{^1) 
is single peaked for E smaller than a critical threshold 
Ec. For E larger than Ec, the distribution P{fl) changes 
to a double-peaked distribution. The analog of the or- 
der parameter ^' which describes the above qualitative 
change in behavior is zero for E < Ec and behaves as 
~ |E - Eel'' for E > Ec, where /3 « 1. We have also 
seen that the dynamics of demand can be understood in 
terms of a 'stochastic' potential which changes its behav- 
ior with E. As E crosses the critical threshold Ec, the 
system undergoes a transition from a 'disordered' state 
with most probable demand equal to zero to an 'ordered' 
state with two phases, excess demand and excess supply. 

We thank the National Science Foundation for support. 
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FIG. 1. (a) Conditional density P(r2|E) for varying E. 
The distribution changes from a single-peaked distribution 
for small E (solid line), to a double-peaked distribution for 
large E (dashed line), (b) Same as (a) in semi- logarithmic 
scale, (c) Order parameter, (position of the maxima of the 
distribution P(n|E)), as a function of E. For small E, P(n|E) 
displays a single maximum whereas for large E two maxima 
appear. Both E and O have been normalized to zero mean 
and unit first moment. For a more accurate estimation of the 
location of the extrema, all densities were computed using the 
density estimator of | 
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FIG. 2. (a) The drift part u(f2) and (b) the variance part 
av{VC) of Eq (^), estimated from the data as the local mean 
and local deviation of Afl. The curves in (b) have been shifted 
vertically for clarity. 
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FIG. 3. Stochastic "force" F{Q,) for three different values 
of E calculated from Eq. ( pj| ) . The derivative a dt; /dfi is 
calculated by first fitting the function crt; by a third order 
polynomial. For E < Ec, the function F{Q,) displays only one 
zero at n = 0. For E ~ Ec, the existing root starts to branch 
into three roots, and for E > Ec, we find three roots: = 
and Q. = Q±. 
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FIG. 4. Stochastic potential V{Q,) for E < Ec shows one 
minimum at = 0, whereas for E > Ec the potential has 
a maximum at n = and two new minima Q± appear. We 
compute Vs{^) by integrating F{Q) [Fig. ^. The curves have 
been shifted vertically for clarity. 
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